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There are two types of edge states in graphene with/without magnetic field. One is a quantum 
Hall edge state, which is topologically protected against small perturbation. The other is a chiral 
zero mode that is localized near the boundary with/without magnetic field. The latter is also 
topological but is guaranteed to be zero energy by the chiral symmetry, which is also responsible for 
massless Dirac like dispersion. Conceptual roles of the edge states are stressed and reviewed from a 
view point of the bulk-edge correspondence and the topological order. 



I. INTRODUCTION 

Many of physical states of matter are characterized 
by order parameters based on the symmetry breaking. 
Especially a symmetry with continuous parameters has 
special importance in a quantum state of matter. When 
the ground state of the system with continuous symme- 
try does not posses the symmetry, that is, the ground 
state is not invariant against the symmetry operation, 
there exists a gapless excitation as the Nambu-Goldston 
mode. A typical example is the Heisenberg magnet with 
antiferromagnetic order ( Neel order ). The ground state 
of the spherical symmetric ( in spin space ) hamiltonian 
is not invariant against the spin rotation. Then above 
the ground state, there is a gapless spin excitation as the 
spin wave. This excitation is realized by infinitesimally 
small spin deformation over infinitely wide range as was 
done by the Lieb-Schultz-Mattis^. 

Although it has been quite successful to characterize 
many of phases, it has been realized that there still exist 
many important physical phases that are not well de- 
scribed by this spontaneous symmetry breaking recipe. 
This class of matter includes many of quantum Hall 
states, Haldane spin chains (integer spin chains), Kondo 
insulators and much more. A ground state of graphene 
also belongs to this new class. In a ground state with 
strong quantum fluctuation, order formation is strongly 
suppressed and a quantum ground state without any fun- 
damental symmetry breaking is realized. In such a case, 
there is no reason to expect gapless excitations, which 
results in a gapped quantum state. Phases of this class 
of matter are the quantum liquids. Then to character- 
ize such a quantum liquid phase is one of the important 
problems. A novel concept as the topological order that 
was first proposed to discuss the quantum Hall stated is 
now under active studies for the purpose. 

Due to be the absence of symmetry breaking, one needs 
to use something new to describe the topological order. 
One of such possibilities to use geometrical phases of 
the quantum states^^^. It works well theoretically to de- 
scribe some of quantum liquids and spin liquids^'^. An- 
other tool is to use edge states, which are one of the fea- 
ture of the topological ordered states. There exist non 
trivial edge states, as the bound states near the bound- 
aries, appear when the bulk are topologically non-trivial. 



Even if the quantum liquid does not posses any sym- 
metry breaking and featureless, there exist characteristic 
boundary states when the system has boundaries and 
impurities. These generic edge states are characteris- 
tic to the topological ordered states and they themselves 
characterize the quantum states, which we call the bulk- 
edge correspondence^. These kinds of non-trivial edge 
states appear in many different ares in condensed matter 
physics, such as optical lattices^ and photonic crystals^. 
Further an entanglement entropy, which is also useful to 
characterize the quantum liquids, is also directly related 
to the edge states^^^^^. In this short review, we have fo- 
cused on non trivial edge states of graphene based on the 
bulk-edge correspondence and try to make clear topo- 
logical aspects of graphene with and without magnetic 
field. 



II. QUANTUM HALL EFFECT OF GRAPHENE 

A. Edge states of the Quantum Hall states 

As is well known, the Hall conductance of two di- 
mensional electrons has an intrinsic topological meaning 
that is an origin of high accuracy of the quantization. 
This quantization is clearly demonstrated by Laughlin 
using an adiabatic process with gauge invariance where 
Aharonov-Bohm(AB) flux $ passes through the hole of 
the cylinder^^ (Fig{T]). Quantized value of the Hall con- 
ductance (in a suitable unit) is a number of electrons 
carried across the system when the AB flux is adiabati- 
cally increased by one flux quantum. Considering a sys- 
tem with edges, this quantized but ambiguous integer is 
uniquely specified by the number of edge modes Ij in the 
j-th Landau gap where the Fermi energy liesi^ji^ 
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It implies that the edge states has essential importance 
for the physics of the quantum Hall effects. 

Although the Hall current is dissipation less, it is nat- 
ural to describe the Hall conductance as a bulk property. 
It is given within the linear response theory a^, cf^y^^ = 

2 

^ Ci where the summation is over the filled Landau 
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FIG. 1: Cylinder used in the Laughlin's argument. Edge 
states in this review are all discussed in this cylindrical geom- 
etry. 
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FIG. 2: (a) One particle energy spectrum of graphene under 
a magnetic field. Horizontal axis is a magnetic flux per 
hexagon in a unit of magnetic flux quantum, (b) Schematic 
Dirac sea without magnetic field and Landau level of the Dirac 
fermions 



levels and Ci is the first Chern number of the ^-th Lan- 
dau band, Q = 2^ / d'^k {rotk Ai)^, = (V^^l V/^V^^), 
where lipi) is a one body wave function of the i-th Lan- 
dau band^^. 

One particle spectrum of the graphene in a magnetic 
field is given in FigO When the Fermi energy is set 
around the E = 0^ one needs to fill the Dirac sea, which 
causes numerical difficulties. In such a case, following non 
Abelian formulation for the Hall conductance is useful 
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where {x^^x^) = (kx^kyY'^'^^ . Especially as for the 
numerical evaluation of the weak field limit. By us- 
ing this method with a technique invented in the lattice 
gauge theory the Hall conductance of the graphene 
as a function of the chemical potential fi is calculated 
without any numerical difficulties^^. By this technical 
but important development, we have further evaluate the 
Hall conductance of realistic graphene with multi-band 
orbital^. The results by the single band case is plotted 
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FIG. 3: Hall conductance of the graphene as a function of the 
chemical potential /i. ( = 1/31) 
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FIG. 4: Energy spectrum of the graphene on a cylinder with 
zigzag edge via momentum along the cylindrical direction ky . 
Red and blue lines are energies of the edge states localized at 
right and left edges. (0 = 1/21) 



in Fig|3]with the density of states (DOS) without mag- 
netic field^^. Three different behaviors ( electron like, 
Dirac like and hole like ) are clearly observed which are 
separated by the singularities of the DOS. Near £^ = 0, 
it reads cr^^^^ = ^{2N + 1) with integer TV. This is a 
Hall conductance of the Dirac fermions. 

Although he Laughlin's argument is enough to guar- 
antee the integral nature of the Hall conductance, only 
the generic consideration can not fix this integer. It is 
specified by considering edge states that is implicitly as- 
sumed in the Laghlin argument The Hall conductance 
is given by the number of edge modes in the energy gap 



that the Fermi energy hes. In FiglH energy spectrum 
of graphene on the cyhnder ( with zigzag boundary ) 
is shown. The red and blue hues in the energy gap re- 
gions are the energies of the edge states locahzed near 
the right and left boundaries. Counting the number of 
edge mode in the energy gap region, we have obtained 
the Hall conductance of the graphene with edges which 
does coincide to that given by the Chern number, eq.([l]). 



B. Bulk - Edge correspondence 

It can be natural to expect that two expressions for 
the Hall conductance by the bulk and by the edges give 
the same value since both describe the same physics. It 
is actually proved even mathematically as -li^ 

Edge ^^Bulk 
^ xy ^ xy 

which is also written as 

This general structure is also applied for the 
graphene^^. As discussed in the previous section, the 
Chern number Cj is a topological quantity. Also the 
number of edge modes, Ij has also topological meaning. 
Here let us give rough idea of the topological meaning of 
the edge mode^. When one discusses usual electrons 
in a magnetic field 5 in a Landau gauge in x-direction, 
two dimensional hamiltonian is decomposed into the sum 
of the one dimensional hamiltonian of the harmonic os- 
cillators (with parameter ky). When one describes the 
Bloch electrons as graphene on the cylinder with the Lan- 
dau gauge, similar decomposition of the two dimensional 
hamiltonian is possible. In this case, corresponding one 
dimensional system with parameter ky is not a simple 
harmonic oscillator but the Harper equationi^ii^. 

H{2D electrons with B) = i^/^^ (harmonic oscillator) 

ky 

i7(graphene with B) = i^/^^ (harper problem) 

ky 

To describe the edge states ( bound states ) and the 
Bloch states ( scattering states ) on the same footing, 
we need to consider a complex energy surface. As for 
the Harper equation, there are multiple (g) energy bands 
that correspond to the Landau levels. We need to pre- 
pare two complex energy planes and R~ by making 
these q energy bands as branch cuts as in Fig|5](a). To 
discuss the problem, it is useful to make them into Rie- 
mann spheres by identifying the infinities into each single 
points. Then the two Riemann spheres are glued into one 
along the branch cuts (Fig|5ja)-(e).) In Figs|5l we have 
shown this procedure symbolically for q = 3 case. Finally 
we have Riemann surface (RS) with g = q — 1 holes ( han- 
dles ), which coincide to the energy gaps of the Harper 
equation. On this RS, generic Bloch function has several 
zeros which correspond to the energies of the edge states. 
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FIG. 5: Construction of the Riemann surface of the Harper 
equation 



One can identify the edge of the cylinder where the edge 
states localize by the position of the zeros. When it is 
on the the edge state is localized at the right 

(left) edge^^. This is for fixed ky. Then by changing ky, 
these zeros of the Bloch state ( edge state energies ) move 
and form loops since ky = and ky = 27r to be identified. 
Then the winding numbers of the loops around the holes 
( gaps ) of the RS give the numbers of edge state, Ij is 
defined. By using the Laughlin's argument, the winding 
number Ij gives the Hall conductance when the chem- 
ical potential lies in the j-the gap^^. In this way, the 
relation eq.([2]) implies a bulk-edge correspondence of the 
topological quantities^. It is clearly demonstrated in the 
graphene quantum Hall effectsi^. 



III. CHIRAL SYMMETRY AND ZERO MODE 
EDGE STATES 

A. Universality of the Dirac fermions and the 
chiral symmetry 

As discussed in the previous section, graphene under a 
magnetic field has edge states, which are topologically 
stable and characteristic to the quantum Hall effects. 
Since the graphene has a chiral symmetry, there are ad- 
ditional edge states, which are protected not only by the 
topological constraint but also by the chiral symmetry. 
The chiral symmetry requires the edge states to have the 
zero energy as discussed below. 

As is well known the energy dispersion of graphene is 
the Dirac fermions like and vanishing linearly at some 
momentum (FiglHl^a)). Another example of the Dirac 
like dispersion in two dimensions is also known in con- 
densed matter. That is the dx2_y2-wdive superconduc- 
tor (Fig|6](b)) which has been discussed for the high- 
Tc materials^^^^^'^^. In momentum representation, the 
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FIG. 6: Typical energy dispersions for (a) graphene and (b) 
two dimensional d^2_y2-waye superconductor. 



hamiltonians for both cases are written as 



H = ^clHkCk 



From this two examples, we have realized that there are 
universality of the Dirac like dispersion when the system 
has a chiral symmetry. This universal feature further 
guarantees the physics of boundary states as discussed 
below. 



B. Berry phases and the Edge states 

The doubling of the generic Dirac fermion with chi- 
ral symmetry is a universal property of the bulk. This 
bulk property also restricts the existence of special edge 
states as an example of the bulk-edge correspondence. 
As for the chiral invariant system, generic condition for 
the existence of the zero mode on a cylinder (FigHJ) are 
give n^^i^^ . When one takes a momentum along the cylin- 
der as ky^ the Berry phase ( Zak phase ) for each ky, 
^(ky)^oi the one dimensional hamiltonian {ky fixed) is 
defined using the Bloch state |V^/e), (Hkl^l^k) = E{k)\il)k))^ 



where Ck is a two component spinor of the fermion op- 



erators. They are respectively given as Ck 



and 



Ck = Hj. = and = for the graphene and 
the dx2_y2-wdive superconductor respectively where 
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As for the graphene, we have assumed that the hopping 
is only allowed between the different sublattices o and •. 
Then the one particle hamiltonian has a chiral invariance 



As for the d^^^ -2/2 -wave superconductivity, when the or- 
der parameter does preserve the time reversal, the order 
parameter, A(/c), is real (A(/c) G M). In this case, it also 
has a chiral invariance 



{if^,r^} = o, r 
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Then the energy dispersion is given by 
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C 3 Z{k) 
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for the graphene 
for the TR invariant d-wave 



It guarantees existence of even number of Dirac fermions, 
since the zero of the Z{k) gives the linearly vanishing 
energy gap (generically) and the closed curve C{ky) = 
{Z{kx, ky)\kx G [0,27r]} moves on the complex Z plane 
under the condition C(0) = C(27r)^^. It implies the curve 
cut the origin in the complex plane even number of times 
when one change ky : ^ 27r. 
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Of course, it does have an ambiguity of the gauge ( phase 
of the Bloch state IV^^), 7{ky) is well defined up to mod- 
ulo 27r.^. Generically speaking the Berry phase can takes 
any values. However the chiral symmetry of the hamilto- 
nian requires the Berry phase has to be quantized into 
or TT. Using this quantized Berry phases of the bulk, one 
can guarantee the existence of the edge states for an in- 
finitely long (along x-direction) cylindrical system when 
the Berry phase jiky) is tt (mod 27r) as far as the chiral 
symmetry is still preserved with edge^. 

j{ky) = TT mod27r^ 
Two zero energy localized modes 

near the at the right and left boundaries of the cylinder 

To carry out concrete evaluation of the Berry phases an- 
alytically, one needs to fix the gauge^. However, this 
gauge fixing is not necessary, at least, to obtain numeri- 
cal value^. 

Since this condition is applied for each momentum ky, 
the zero energy edge states form doubly degenerate flat 
bands when one considers the original two dimensional 
problem on the cylinder (FigH]). Generically speaking, 
this condition is not satisfied for all ky, then the flat 
bands only exist for restricted momenta. Applying this 
generic consideration for graphene, one predicts flat band 
zero energy edge states for 1/3 of the total momentum 
along y direction near the zigzag edge (FigCj)^^, which 
was found by fujita et al.^^. This edge mode does not 
exist for armchair edges. It is also consistent with this 
generic condition (Fig|7j). 

As for the dx2-y2-wdive superconductor, we have ap- 
plied it and obtain the flat band zero mode edge states 
for (llO)-direction^^. As for the (lOO)-direction, we pre- 
dict absence of the zero energy edge states^^ ( There is 
no topological reason to have the edge states. See FigCj). 
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FIG. 8: Energy spectrum of graphene under magnetic field 
on the cy Under with zigzag edges. The flux per hexagon is 
(j) = 1/4, 1/7 and 1/21 in unit of flux quantum respectively. 
Red and blue lines are the edge states localized near the left 
and right edges respectively. The gray regions are the Landau 
levels. 



FIG. 7: Energy spectrum of the graphene and (ia,2_^2-wave 
superconductors with different angles on cylinders. 

This is consistent with the existence of the Andreev 
bound states only for the (llO)-direction, which has been 
observed in the zero bias conductance measurement^. 



C. Spontaneous local chiral symmetry breaking 

As discussed, graphene and the dj;2_y2-wdive supercon- 
ductors have zero energy edge states protected by the 
chiral symmetry. It is topological in the sense that the 
quantized Berry phase ( Zak phase ) is a key (bulk) quan- 
tities to guarantee the edge states. 

These zero energy edge bands are completely flat ( 
^ = ) if they exist. The flatness is guaranteed by 
the chiral symmetry. Of course, the symmetry of the 
bulk can not be changed by the existence of the bound- 
aries. However, edge potential induced by the edges may 
break the bulk symmetry. This symmetry breaking is 
allowed to exist only near the boundaries (by the edge 
potential). When one regards the flat band zero modes 
as a one-dimensional system, it is natural to expect the 
Peierls instability that reduces the symmetry since the 
flat bands have a diverging density of states. Then what 
would be expected is that the chiral symmetry has bro- 
ken only near the edges to make the flat band dispersive 
or making a gap in the edge mode bands. This spon- 
taneous symmetry reduction occurs only local near the 
boundaries'^. This actually occurs both in the graphene 
and the d^^^-ys-wave superconductors. In the case of 
graphene, it corresponds to the boundary magnetic mo- 
ments near the zigzag edges. When one considers this 
boundary magnetism within the mean field theory, local 
moments induce site dependent potentials. It destroys 
the chiral symmetry of the graphene near the bound- 



aries. It is confirmed, at least, numerically based on the 
density functional theory calculation'^. 

The same local symmetry breaking also occurs in the 
dx2_y2-wdive superconductors^^. In this case, the chiral 
operator corresponds to the time-reversal operation that 
implies real order parameters of the superconductors. 
Then what happens is that the order parameter of the 
superconductivity becomes complex only near the bound- 
aries. It implies spontaneous generation of local fluxes. 
Topologically, the boundaries of the cylinder and point 
like impurities are equivalent. Therefore, the generic con- 
sideration predict that spontaneous local flux generation 
near the boundaries or impuritie o^^i^' . 



D. Chiral zero modes under the magnetic field 

When one considers the graphene under magnetic field, 
there are two different type of edge states. The one is 
topologically protected quantum Hall edge states. The 
other is chiral zero mode edge states, which are topo- 
logical, but are protected by the chiral symmetry. The 
former is discussed in the previous section. Then we dis- 
cuss the latter one here. 

As shown in Fig|8l there co-exist Dirac fermions' n = 
Landau level and the chiral zero mode edge states for 
graphene with zigzag edges. Without magnetic field, 
there also exist chiral edge modes. However, it is de- 
generate with the bulk Dirac fermion at the gap closing 
momentum. Therefore there is no length scale for the 
E = states. With the magnetic field, the situation is 
different and the length scale should be determined by 
the magnetic length scale. This additional contribution 
of the boundary state to the local charge occurs only 
near the zigzag edge. It is absent at the armchair edges 
(See Figl9])^^. Generically speaking, the edge potential 
repels the charge density to make depletion layer near 
the boundary Physically the length scale of the deple- 
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FIG. 9: Local density of states of the graphene with several 
edges around ^ = 0. Distance from the edge is denoted by x 
and is a magnetic length. (a):zigzag edge and (b) armchair 
edge. 



out magnetic field. However with magnetic field, there 
also exists a n = Landau level of the Dirac fermion at 
E = 0. It plays fundamental role in reconstruction of 
the boundary charge. Local charge density of this ^ = 
states, which has two physical origins as the chiral edge 
states and the n = Landau level, has been investi- 
gated by Arikawa et al^'*^ They have calculated the lo- 
cal density of state. It can be experimentally observable 
by the scanning tunneling microscope (STM). Although 
there exist chiral zero modes at the zigzag edge, it is 
not enough to describe the edge charge reconstruction. 
The chiral zero modes actually show detailed structure 
and long tail into the bulk. However, what is realized in 
the total local charge density is exponentially localised 
structure and it becomes to the bulk value quite rapidly. 
It will be understood as the topological compensation 
(screening) of the chiral edge charge by the bulk. This 
is a new feature of the boundary reconstruction of the 
^ = states of graphene. 



tion layer is determined by the electro-static potential of 
the edge detail. Without including it, the length scale 
is the magnetic one is- It is realized at the armchair 
edge.(See[9]:b)) 

The situation is different at the zigzag edge, there oc- 
curs enhancement of the local charge density instead of 
the depletion as shown in Fi^a). This is consistent with 
the existence of chiral zero modes, which also exist with- 
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